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We present methods which allow orders of magnitude increase in the number of modes in linear 
optics experiments by moving from spatial encoding to temporal encoding and using dispersion. This 
enables significant practical advantages for linear quantum optics and Boson Sampling experiments. 

Passing consecutively heralded photons through time-independent dispersion and measuring the 
output time of the photons is equivalent to a Boson Sampling experiment for which no efficient 
classical algorithm is reported, to our knowledge. With time-dependent dispersion, it is possible 
to implement arbitrary single-particle unitaries. Given the relatively simple requirements of these 
schemes, they provide a path to realizing much larger linear quantum optics experiments including 
post-classical Boson Sampling machines. 


Unitary transformations on optical modes have been 
used to implement single particle quantum gates [1, 2], 
quantum simulations [3], and Boson Sampling [4-11]. 
Traditionally, these transformations are implemented on 
spatial modes using a system of beamsplitters. How¬ 
ever, building a large interferometer implementing such a 
unitary transformation is experimentally challenging and 
the largest number of modes so far has been 21 [9]. 

A particularly interesting application which requires a 
large number of spatial modes is Boson Sampling. Boson 
Sampling is the process of estimating the output pho¬ 
ton distribution after passing multiple identical photons 
through a passive linear interferometer. Aaronson and 
Arkhipov have proposed that Boson Sampling is compu¬ 
tationally hard for classical computers because it requires 
the estimation of the permanents of independent and 
identically distributed (iid) Gaussian matrices, a prob¬ 
lem that is believed to reside in the #P-complete com¬ 
plexity class [4]. The only ‘effective interaction’ between 
photons in such a system is due to the bosonic statis¬ 
tics of these identical photons at the detectors. The need 
for only linear optics and photodetection could make the 
Boson Sampling problem easier than general quantum 
computing approaches with photons, which require non¬ 
linear materials [12] or feed-forward schemes [13]. Fur¬ 
thermore, unlike other quantum computing schemes that 
require on-demand sources, Boson Sampling with proba¬ 
bilistic but heralded input photons has been proposed to 
be computationally hard for a classical computer [14]. 

In conventional Boson Sampling schemes that use spa¬ 
tial modes (which we now refer to as ‘Spatial Mode Boson 
Sampling’ or SMBS), multiple identical photons enter 
a high-dimensional transformation over spatial modes, 
such as a system of beamsplitters and phase shifters, 
while the output probability distribution is monitored 
with detectors at each of the output modes [5-11], as 
shown in Fig. 1 (a). Specifically, photons are injected 
into input modes 1 The system then trans¬ 

forms the creation operator for input spatial mode j as 
CLj -T JA Uij&l, at which point photons in each mode are 


measured using single photon detectors [4-9]. As we de¬ 
scribe, Boson Sampling can analogously be performed in 
time by replacing spatial mode aj with temporal mode 
a t j (Fig. 1 (b)) and by replacing the beamsplitter array 
with dispersion (Fig. 2). 
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FIG. 1. Schematic of (a) spatial mode Boson Sampling 
(SMBS) (b) temporal mode Boson Sampling (TMBS) 

SMBS entails several difficult challenges, that, as we 
show, favor temporal mode encoding. First, Boson Sam¬ 
pling requires an extremely large number of modes to be 
classically computationally difficult. Strictly speaking, 
the complexity argument for Boson Sampling assumes 
that, if n is the number of photons in the system and 
m is the number of modes, m > Q(n 5 log 2 n). Although 
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FIG. 2. Schematic of TMBS with time-independent disper¬ 
sion 

Aaronson and Arkhipov have conjectured that the com¬ 
plexity arguments still hold when m = 0(n 2 ) [4], the 
number of modes is still large: e.g., even with m = n 2 
and n = 30, the interferometer would require 900 modes. 
Furthermore, the experiment would require 900 detec¬ 
tors and, if the photons came from heralded sources, 900 
sources [14]. To date, experimental demonstrations of 
SMBS have been limited to 5 photons in 21 modes [9]. 
The number of modes in temporal mode Boson Sampling 
(TMBS) can be increased simply by increasing the dis¬ 
persion. Even with 10000 ps/nm dispersion, which can 
be achieved with off-the shelf components, and 100 ps de¬ 
tector jitter, which can be routinely achieved with silicon 
avalanche photodiodes or with superconducting nanowire 
single photon detectors [15], the number of modes in 
TMBS is orders of magnitude higher than SMBS. In prin¬ 
ciple, the experiment can be implemented in a single fiber 
with only a single photon source and only two detectors: 
one to herald input photons and one to detect the output 
state, regardless of the number of interfering photons in 
the system. Given the dead time tdt of single photon de¬ 
tectors, the output may have to be split between a larger 
number of detectors. However, in general, the number of 
detectors is smaller than required in SMBS [16]. 

Furthermore, uncertainty in the time when photons 
are injected into different modes leads to distinguisha- 
bility and loss of boson interference. This is a particular 
problem in SMBS with heralded sources based on sponta¬ 
neous parametric down conversion (SPDC). Most SMBS 
experiments to date have relied on downconverted pho¬ 
tons. Temporal or spectral filtering could improve the 
interference, but at an exponential loss in multi-photon 
throughput. In TMBS, the lack of control over the input 
time of our photons only corresponds to a lack of con¬ 
trol over the choice of our input modes; however, as long 
as the input modes are known, this does not affect the 
ability to perform Boson Sampling [14]. 

Previous proposals have considered temporal modes 
for Boson Sampling [17, 18] but they relied on temporar¬ 


ily converting temporal modes to spatial modes and 
then mixing the modes with beamsplitter operations. 
Hence, increasing the number of output modes (m n) 
requires a large number of effective beamsplitter oper¬ 
ations. They also require active elements that operate 
on a picosecond time scale. Furthermore, since they are 
based on the interference of narrow photon packets, they 
suffer from the same issues with temporal mismatch as 
SMBS. In TMBS, detector jitter can limit the accu¬ 
racy with which the input mode can be heralded but 
this limitation can be overcome by using large dispersion. 

Time-independent dispersion - We first consider the sim¬ 
plest case of Boson Sampling in time with identically 
shaped input photons and time-independent dispersion 
(Fig. 2). If we use an SPDC source with idler photons 
heralded at times tj and signal photons used as input 

photons, the input state is given by |\Eq n ) = (n, A) i°) 
where |0) is the multimode vacuum state and a\- = 
fZdt a)(t)A(t — tj) represents the creation operator 
for the input state centered at tj. cA(t) is the creation 
operator for time t and cjo is the central frequency of the 
input photons. We assume that the photon state after 
heralding of the idler is a pure state of the form A{t — tj). 
However, a realistic detector projects the signal photon 
into a mixed state with tj varying over the timescale of 
the detector jitter; the effect of this temporal mismatch 
can be made negligible with large dispersion [16]. 

a\- can be expanded in the frequency domain as 
d\- = dujcA {uj)^ {A{t — tj)} where at(u;) is the cre¬ 

ation operator for frequency uj and iF{A(t — tj)} is the 
Fourier transform of A(t — tj). After passing through 
a dispersive element with dispersion relation /3(uj) and 
length L, frequency components at u) are multiplied by 
a factor where <fi(uj) & f3(uj)L. The wavefunction 

of the multi-photon system is then given by | 4/ out ) = 

(llj |0> where fot = dujd\uj)J r {A{t — tj)}e~ l ^ UJ \ 

Going back to the time domain, V- = d tcA {t)U (£, tj) 
with 

U(t,t j )=A(t-t j )*F- 1 {e ~ i *">} (1) 

where V is the convolution operator. 

If dispersion parameters are chosen such that U(t,tj) 
does not change appreciably when t varies in a window 
of width t s [16], the modes can be discretized so that the 
transformation is well approximated by a\- —>• JT Uijd ^ 
where cA ti represents the creation operator at the dis¬ 
cretized time step near G, i.e. a\- = J^ ts dtaA (t) / y/t^. 
We can then write the transformation as a\- —>> 

JT Uijd^ = dtU(t : tj)a^(ti). If we assume that 
U(t,tj) is approximately constant for a small time step 
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£ s , we can approximate Uijd) ti ~ U(ti,tj ) f* J+ts dta^(t). 

Hence, we have Uij = y/t^U(ti,tj). d\- 
where 


Uii = 


VT s A{t-t 3 :)*^- 1 {e" i ^“)} 


( 2 ) 


Eq. 2 shows the class of unitary transformations from 
which we can sample using time-independent disper¬ 
sion. The unitary is band-diagonal because of the time- 
invariant nature of the system. Classical algorithms exist 
for the computation of the permanent of banded matri¬ 
ces with a banded inverse which is polynomial in the size 
of the matrix but exponential in the number of bands 
[19]. The inverse of a unitary banded matrix is banded. 
However, because the number of bands is extremely large 
and the number of bands/dispersion is increased with the 
number of photons in order to limit the effect of jitter 
[16], the problem is still expected to be computationally 
hard. 

The shape of the input pulses A{t) is incorporated into 
the unitaries in Eq. 2 because, unlike conventional uni¬ 
tary implementations, the input states and measurement 
have different bases; the input photons have shape A(t) 
but the measurement is in the time basis (with eigen¬ 
functions S(t — tj)). The results of an experiment will 
be the same as a spatial unitary implementing Eq. 2. 
Imagine a fictitious experiment where the input photons 
are 5(t — tj). They then go through a unitary U\ which 
puts them in a superposition of the form A(t) (physically, 
this is the state of the photons going into the dispersion). 
The photons then go through another unitary U 2 which is 
the dispersion and the total unitary implemented by this 
system is UiU 2 - The output of this system will be the 
same as our scheme. Although a wavefunction S(t — tj) 
is unphysical, the detector sees the same output state as 
if the operator U 1 U 2 was applied to photons of the form 
5(t — tj). 

It is possible to sample from a larger class of uni¬ 
taries by shaping the temporal form of the input pho¬ 
tons. Methods for shaping single photons with arbitrary 
amplitude and phase in time have been proposed [20] 
and an experimental demonstration of shaping the spa¬ 
tial waveform of single photons had been reported [21]. 
With pulse-shaping, the input waveform A(t — tj) is re¬ 
placed by a more general set of functions Aj(t) so that 
the accessible set of unitaries becomes 


Uij 




t=ti 


( 3 ) 


If it were possible to choose any set of functions Aj (t ), 
then the unitary could be chosen column by column using 
Aj and simply detecting the photons without any disper¬ 
sion would be equivalent to Boson Sampling. However, 


it is experimentally challenging to prepare multiple over¬ 
lapping photons with a specific waveform. Hence, for re¬ 
alistic implementation, the photon wave-packets should 
be separated in time. This limits the possible unitaries 
represented by Eq. 3. 

Although we have no proof of the hardness of sampling 
from such a unitary, there is to our knowledge no reported 
efficient classical algorithm for sampling from a general 
unitary of this form. Sending a train of photons through 
a time-independent dispersion could allow for a Boson 
Sampling experiment with more photons and more modes 
than can be currently achieved with SMBC. 

An arbitrary functional form for the dispersion 0(u;) 
can be obtained by using approaches used in optical func¬ 
tional design [22] and femtosecond pulse-shaping [23]. 
There are even commercial products for implementing 
arbitrary dispersion used for pulse-shaping in telecom¬ 
munication [24]. 

A central feature in spatial boson collision experiments 
is the bunching of bosons in the output modes i.e., Hong- 
Ou-Mandel interference [25]. An analogous feature ap¬ 
pears in the temporal modes. Consider a heralded in¬ 
put state of two photons, (A (t — t\) a) (t — t 2 ) |0). If this 
state passes through group velocity dispersion in a fiber 
of length L : 0(c j) = j8 0 + j8i(cj — cj 0 ) + 1/2 f3 2 (uj - cj 0 ) 2 , 
4>j = PjL,j G {0,1,2}. For simplicity, we have assumed 
a Gaussian input shape of the form 


m = 



1/2 

exp 


■_ t^J 

4<t^ 

L cor J 


exp[iu)ot\, 


( 4 ) 


where we have assumed that the correlation time a cor 
of the photons from the heralded photon source is much 
shorter than the biphoton coherence time. SPDC pho¬ 
tons are often approximated as Gaussians in time. How¬ 
ever, the two photon interference effects would be visible 
with any shape in general. After passing through this 
dispersive element, the probability of detecting the pho¬ 
tons at times t\ and t 2 corresponds to the magnitude 
squared of a 2 x 2 permanent and assuming 0 2 cr? or , 
the probability goes to zero when [16] 


At l At° = (2 m + 1)7t02 (5) 

In Fig. 3, the joint probability of observing a photon 
at t\ and t 2 is plotted when two photons with cr cor = 
200 fs and a Gaussian temporal waveform centered at 
t\ = —0i and t 2 = — 0i + 100 ps are sent through a 
dispersive element with a GVD parameter of magnitude 
\D\ = 2iYC<fi 2 /\ 2 = 10000 ps/nm. The assumption 
02 tf'cor h as n °t been used. At the time-scale of a 
few nanoseconds, the output resembles a Gaussian pulse 
centered at t° = 0, as would be expected from single 
photon input at t l = —0i- However, at the time scale of 
hundreds of picoseconds, two-photon interference effects 
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FIG. 3. (a) Setup for seeing ’HOM-like’ interference (b) The 
joint probability of detecting the first photon at t° and the 
second photon at when two input photons near t = 0 and 
separated by 100 ps are sent through a second order disper¬ 
sive element. cr cor = 200 fs and the dispersive element has 
a GVD parameter of magnitude \D\ — 27TC02/A 2 = 10000 
ps/nm. t s — 10 ps 

can be seen in clear dips in the two-photon output 
probability appear, as predicted by Eq. 5. The plot has 
been discretized with t s = 10 ps. We chose a low value 
of t s here to clearly show the shape of our interference 
pattern; such timing resolution is not necessarily re¬ 
quired to resolve the two-photon interference pattern. 
Experimentally, these dips would be easily resolved 
using detectors with a jitter of 100 ps. The predicted 
correlation measurement is provided in Fig. 1 of the 
supplemental material [16]. The increase in the size of 
each bin would also increase the probability of detecting 
photons in each time bin by a factor of 100. 

Time-dependent dispersion - With time-dependent dis¬ 
persion, it is possible to implement any arbitrary unitary 
transformation on temporal modes. 

Using the time projection operator t p = / d t't \t') (t'\ 
and the frequency projection operator ti p = 


FIG. 4. Schematic of TMBS with time-dependent dispersion 

fdcj (j\uj)(u\ which satisfy [t p ,ti p \ = i/2 (similar 

to x and p in [26, 27]), an arbitrary continuous variable 
unitary on temporal modes can be written as e~ l ^ tp,UJp ^ 
where f(t p ,ti p ) = f dtH(t p ,ti p ,t). H{t p ,ti p ,t) is the 
Hamiltonian of the system. It should be noted that 
t p is the time projection operator corresponding to the 
photon’s temporal waveform and not the time over which 
the system evolves. Hence, the variable of integration t 
above is distinct from the operator t p . 

Such a Hamiltonian can be physically realized by mak¬ 
ing the elements of the dispersion time-dependent. In the 
previous section we discussed methods for implementing 
arbitrary dispersion which implements a Hamiltonian of 
the form 4>(ti p ). Since photons reaching the dispersive el¬ 
ement at different times see a different 4>(ti p ), the Hamil¬ 
tonian can be written as H(t p ,ti p ,t) where H can be 
any real function. Hence, we get the desired unitary 
e -tf(tp,u p ) w ith time-dependent elements that allow for 
changes in the frequency spectrum of the pulses which 
was not possible with time-independent dispersion. 

Another option for realizing such a unitary would be 
to cascade elements that implement time-independent 
dispersion and a time-dependent refractive index. The 
Hamiltonians corresponding to dispersion and time de¬ 
pendent refractive index are 4>(ti p ) and g(t p ) respectively. 

Following previous work on realizing arbitrary Hamil¬ 
tonians for continuous variable systems [26, 27], one can 
construct Hamiltonians of the form [gi(t p ),(/)i(ti p )] + 
[ 92 (t p ), [g 3 (ip),<t> 3 (iOp)]} + . using the properties 

e -iASt e -iBSt e iASt e iBSt _ e [A,B]St 2 + 3) (g) 

e iA8t/2 e iB8t/2 e iB8t/2 e iA8t/2 _ Q i{A+B)8t (J(5t 3 ). (7) 

The following Hamiltonians can be cascaded to gen¬ 
erate a Hamiltonian that can be any polynomial in t p 
and Cj p [26]: t p , Cj p , - j- Cj p and a Hamiltonian of the 

form where n > 3. The Hamiltonians ti p and t p are 
first order dispersion (inverse group velocity) and a linear 
varying refractive index. A high order Hamiltonian ti p 
can be realized with higher order dispersion. t p + ti p can 
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be realized with a combination of second order disper¬ 
sion and a quadratically varying refractive index using 
Eq. 7. Although such a decomposition allows us to build 
arbitrary Hamiltonians with a number of elements which 
increase as a small polynomial in the number of photons 
[28], more efficient decompositions are often possible with 
fewer elements [27]. 

Furthermore, as opposed to the conventional construc¬ 
tion of f(x,p) [26] which uses low-order polynomials in 
x and p (high order polynomials require a non-linear 
medium), a unitary of the form g{t p ) or (j){Qj p ) can be 
of any arbitrary functional form without requiring an ex¬ 
plicit Kerr-type nonlinearity. 

Given the ability to realize arbitrary continuous vari¬ 
able single particle unitaries over t p , any discrete unitary 
transformation can be implemented by making the trans¬ 
formation constant over the output time bins [16]. Such 
an experiment with multiple photons would be equivalent 
to Boson Sampling and if the discrete unitary is chosen 
with Haar measure, the results are believed to be classi¬ 
cally intractable [4]. 

In conclusion, we have introduced new methods of im¬ 
plementing unitary transformations on temporal modes 
based on dispersion and pulse shaping that require a 
much smaller number of sources and detectors and do 
not require a large system of beamsplitters. In principle, 
using only fixed dispersion, a single heralded source and 
two detectors, one can observe multi-photon interference 
and perform a Boson Sampling experiment for which no 
efficient classical algorithm is known, to our knowledge. 
By using time-dependent dispersion, it is possible to sam¬ 
ple from arbitrary unitaries. 
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Supplemental Material: Scalable Boson Sampling Schemes by interference in the 
Time and Frequency Bases using dispersion and pulse shaping 


NUMBER OF DETECTORS 

In this section we show that the number of detectors 
in TMBS can be much smaller than SMBS 

In Fig. 3 in the main paper, if photon detection is 
binned in steps of 100 ps, there are 1934 modes. The 
number of modes is defined as the number of time bins 
within which the absolute value of the dispersed wave- 
function is greater than 90% of the peak value. 

We assume that our detectors have a dead time of 1 
ns and look at the failure rate of our boson sampling 
scheme with 2000 input and output modes and 30 pho¬ 
tons when each photon output is passively and equally 
split between 30 detectors. We assume that each time bin 
on the heralding as well on the output detector bank is 
equally likely to receive a photon. We post-select on the 
cases where there are a total of 30 photons incident on 
each detector bank (as in conventional scattershot boson 
sampling). 

The scheme is considered a failure if two photons are 
incident on any detector within the dead time. Based on 
the Monte-Carlo simulation of the system, we find that 
the probability of failure is less than 10% 

The assumption of photons being equally likely to ar¬ 
rive at any bin is accurate for the heralding detector 
bank. For the detector bank which detects photons after 
going through the unitary, the assumption may lead to 
an underestimated failure rate. However, we can postse¬ 
lect on the number of photons detected after the unitary 
being equal to the number of heralded photons and hence 
a higher failure rate is tolerable. An accurate simulation 
of the failure rate with a 30 photon, 2000 mode system 
is expected to be close to the limit of current computing 
capabilities. 

Hence, TMBS can allow for a 30 photon 2000 mode ex¬ 
periment with 60 detectors, whereas an equivalent SMBS 
experiment would require 4000 detectors. It is interesting 
to observe that for the same number of photons and dead 
time bins, the number of detectors required for TMBS 
goes down with an increase in the number of modes since 
there is a smaller chance of detecting multiple photons 
within a dead time reduces. In SMBS, the number of 
detectors is equal to twice the number of modes. 

ERROR BOUNDS 

We find bounds on the error in the sampling distribu¬ 
tion due to detector jitter and discretization. U is the 
ideal unitary that we wish to implement, U is the uni¬ 
tary with errors and Vjj and Dp are the corresponding 


probability distributions over outcomes. 

It has been shown in [1] that if there are n photons in 
the system, 

WV^-VuWKnWU-UWop ( 1 ) 

If the relative error of each matrix element has a upper 
bound of R i.e. | Uij — Uij\ < R\Uij\ for all i, j, 


WVo-VuW^nWRUWop 
\\V 0 -Vu\\<nR (2) 

where we have used the fact that U is unitary. 

Hence, in order to have ||Dp — Vjj || = o( 1), R = o(^). 

Error due to detector jitter 

We have shown that the in the case of input photons 
with a fixed shape A(t — tj ) sent through a dispersion 
P(uj)L, the resulting unitary sampled from is 

Uij = (3) 

= y/t a A{U - tj) * B(ti) 

= y/UA(U) * B(ti - tj) 

where B(t) = •7 r-1 {e - * < ^ w )}. 

Due to detector jitter, there is an uncertainty in tj. If 
the maximum timing error due to detector jitter is £ e , 


R = 

A(U) * B(ti - tj + t e ) - A(ti) * B(ti - tj) 

A(ti ) * 

B{ti tj) 


A(ti)*B(ti -tj) 
A(ti) * B(U - tj) 





If we write the dispersion in the form 

4>{u) = <t>'(<t> a u) (5) 

where the (j) s is used to scale dispersion. It can be seen 
that B(ti — tj)/B(ti — tj) = o(l/(f) s ). Hence, if A{t) is 
chosen independent of n, R = o(t e /(f) s ). Therefore, in 
order to limit the error due to detector jitter, t e /<f) s = 
o(l/n). 
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Error due to discretization 

In the case of time independent dispersion, the wave- 
function generated after dispersion is treated discretely 
in order to draw a parallel with the original formalism 
for Boson Sampling [2]. 

Uij = yJtsU(ti,tj). The relative error can be written 
as 


R < max t 


: max t 



(6) 


where £ G [U — t s /2,ti + t s / 2]. If we define ( = 
ma x t {U(t,tj)/U(ti,tj)}, we can see that R = o((t s ). As 
discussed previously, we have ( = o(l/0 s ). Hence, in 
order to limit the error, t s /<f) s = o(l/n). 


TEMPORAL ANALOG OF THE HOM DIP 


In this section, we derive the two-photon wavefunction 
obtained on passing two photons with Gaussian envelopes 
through second order dispersion which is used to derive 
Eq. 5 in the paper. 

Using Eq. 1 of the main paper, a sin¬ 
gle photon with temporal waveform A(t) = 
l 1/2 r + 2 i 

exp[icj 0 £] on passing through 


( i y 

\ & cor y/ 27T / 


exp 


'4 0-2 


second order dispersion results in the temporal waveform 


m = (2 tt)- i/4 


exp 


L 4 (^ 


°cor + *<fe/2 

ft ~ (f>l ) 2 

■4 

cor 


e -i4>o e iuot ( 7 ) 


>1/4) 


.02 


Hence, for two sets of entangled photons generated 
from a heralded source with idlers are detected at times 
t\ and £ 2 , after passing through the dispersive element, 
the wavefunction is given by 


/ oo poo 

d t' / d t" a\t r )a\t") 

-00 J —00 

A(t' -4) (8) 

Every combination of g0(£i)g 0(£2) is repeated twice 
under the integrals. The repetition can be removed by 
rewriting the expression as 


/ OO POO 

d t' / d t"a) {t')a) (t") 
-00 J t' 

| n , _+ *^2/2 


<re 


t iuj 0 (t' +t") — z20 o 


1 + exp 


4(7* 


v^(^ or + i<t> 2 / 2 ) 

2 {2fti-4)ft'-0}j} 


exp 


~Vcor + *^2/2 


^ 4<j4 or + 02 

+ ft"-4-^) 2 }l |o) 


{ft' -4 -<M 2 


(9) 


assuming that the pulse broadening due to the dis¬ 
persion is much greater than the correlation time of the 
photons from the entangled source (0 2 cr;? or ), this is 
reduced to 


/ oo poo 

d t' / d t" a) (t r )a) (t rr ) 

-00 J t' 


|l + exp 

o 

[ 


M-^w- n 

h 


]} 


exp 


4cr cor+ _ 

+ ft" ~ 4 - 0i) 2 }] |0) 




+t")-i2(f)o 

V^^ cor + i<t> 2/2) 


2 


(10) 


From the equation above, Eq. 5 from the paper follows 
immediately. 

Using a], = J^ 3+ts d to){t)/y/t~ s to discretize Eq. 9, we 
get 


l*) = E E «UUer(M) |0) (11) 

tl t 2 >=ti 


where 


Mi 


£ s cr 


1/2 


jk 


p^O tj —<f>o 


Vv^F(^+i0 2 /2), 


exp 


—< 7 ; 


2 

cor 

2 

cor 


( 12 ) 


+ £02/2 


4(7^ 

^ cor 


(ft - 4 


Hence, the probability of detecting photons at £° and £2 
is given by [|Per(M)| 2 /(r in \r out \)] t =t o where r in (r out ) = 
2 if the the input (output) photons are in the same mode 
and 1 otherwise. 

The joint probability of detecting two photons at £° 
and £2 is plotted in Fig. 3 in the paper with t\ = — 0i, 
t\ = —0i + 100 ps, \D\ = 27TC02/A 2 = 10000 ps/nm and 
a cor = 200 fs. In the main paper the detection time has 
been binned in 10 ps steps which is hard to achieve be¬ 
cause of detector jitter. In Fig. 1 here, the two photon 
interference pattern is visible even when the binning is 
increased to 100 ps which is much easier to achieve ex¬ 
perimentally. 
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FIG. 1. The joint probability of detecting the first photon 
at t° and the second photon at when two input photons 
near t = 0 and separated by 100 ps are sent through a second 
order dispersive element. a CO r = 200 fs and the dispersive el¬ 
ement has a GVD parameter of magnitude \D\ = 27TC02/A 2 = 
10000 ps/nm. The probability has been binned into buckets 
of 100 ps which corresponds to a jitter achievable with cur¬ 
rently available silicon and superconducting nanowire single 
photon detectors [3] 





